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Abstract 

We study the perturbations of two classes of static black ellipsoid solu- 
tions of four dimensional vacuum Einstein equations. Such solutions are de- 
ep ' scribed by generic off-diagonal metrics which are generated by anholonomic 

transforms of diagonal metrics. The analysis is performed in the approximation 
of small eccentricity deformations of the Schwarzschild solution. We conclude 

■ that such anisotropic black hole objects may be stable with respect to the per- 
turbations parametrized by the Schrodinger equations in the framework of the 

fSJ ' one-dimensional inverse scattering theory. 

o ■ 
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■ 1 Introduction 

I In a series of works new classes of solutions of four dimensional (4D) and 5D vacuum 

Einstein equations with ellipsoid and toroidal symmetry are constructed P . Such solu- 
tions are generated by anholonomic deformations of the Schwarzschild metric, of static 
or stationary configurations, and depend anisotropically on angular coordinates. They 
are described by generic off-diagonal metric ansatz which are effectively diagonalized 
with respect to anholonomic frames with associated nonlinear connection structure. A 
study of horizons and geodesic behaviour ^ concluded that for small deformations 
of the spherical symmetry, for instance, to a resolution ellipsoid one. Such solutions de- 
fine black ellipsoid objects (static black holes with ellipsoidal horizons and anisotropic 
polarizations of constants). With respect to anholonomic frames the new solutions 
are given by certain metric coefficients which are similar to those from the Reissner- 
Nordstrom metric, but with an effective, polarized, "electromagnetic" charge, induced 
by off-diagonal vacuum gravitational interactions. This differs substantially from the 
usual static electrovacuum solutions with spherical symmetry which are exact solutions 
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of the Einstein-Maxwell equations (see, for instance Ref . j3] , as a general reference on 
"the mathematical theory of black holes"). 

The aim of this paper is to study perturbations of black ellipsoids and to prove 
that there are such static ellipsoid like configurations which are stable with respect 
to perturbations of a fixed type of anisotropy (i. e. for certain imposed anholonomic 
constraints). The main idea of a such proof is to consider small (ellipsoidal, or another 
type) deformations of the Schwarzschild metric and than to apply the already developed 
methods of the theory of perturbations of classical black hole solutions, with a re- 
definition of the formalism for adapted anholonomic frames. 

The theory of perturbations of the Schwarzschild spacetime black holes was initiated 
in Ref. 0, developed in a series of works, e. g. Refs [HIIZI, and related [S] to the theory 
of inverse scattering and its ramifications (see, for instanse, Refs. [0]). The results 
on the theory of perturbations and stability of the Schwarzschild, Reissner-Nordstrom 
and Kerr solutions are summarized in a monograph |3]. As alternative treatments of 
the stability of black holes we cite in Ref. 

The paper has the following structure: In Sec. 2 we introduce an off-diagonal ansatz 
which parametrizes two classes of anholonomic static deformations of the Schwarzschild 
solutions which describe black ellipsoid like objects (the formulas for the components of 
Ricci and Einstein tensors are outlined in the Appendix). In Sec. 3 we investigate the 
axial metric perturbations governed by some one-dimensional Schrodinger equations 
with nonlinear potential and anisotropic gravitational polarizations. Section 4 is de- 
voted to a stability analysis of polar metric perturbations; a procedure of definition of 
formal solutions for polar perturbations is formulated. In Sec. 5 we prove the stability 
of static anholonomically deformed solutions of the Scwarschild metric with respect to 
perturbations treated in the framework of the inverse scattering theory based on the 
one-dimensional Schrodinger equation with paremetric potentials. A discussion and 
conclusions are contained in Sec. 6. 

2 Metrics Describing Perturbations of Anisotropic 
Black Holes 

We consider a four dimensional pseudo-Riemannian quadratic linear element 



'2tti s ' ' • 

1 1 — - ) dr^ — r^q{r)d6^ — rj^{r, 9, ip)r^ sin^ 96ip^ 



+ 

with 



2m e , , 
1 + ^r]{r,Lp) 



6t' 



5ip = dip + ewi{r, (p)dr, and 6t = dt + eni{r, (p)dr, 

where the local coordinates are denoted u = = {r,9,ip,t)} (the Greek indices 
a, 13,... will run the values 1,2,3,4), e is a small parameter satisfying the conditions 
< e <^ 1 (for instance, an eccentricity for some ellipsoid deformations of the spher- 
ical symmetry) and the functions Q{r,ip),q{r), rj^{r,9,ip) and rj{9,ip) are of necessary 
smooth class. The metric (Q) is static, off-diagonal and transforms into the usual 
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Schwarzschild solution if e — and g, 773 1; it describes at least two classes of static 
black hole solutions generated as small anhlonomic deformations of the Schwarzschild 
solution IIII2II31- 

The geodesic and the horizon structure of the first class of solutions was investigated 
in Ref. P). It was proved that for the data 

?74(r,v9) = 1 + e\4{r,(p) - e^\4{r,(p) il - —j + 
r]{r,(f) = X4{r,ip) (r'^ — 2mr^ + 1,Q = 1, 



r]od^/\ru\/d(p, t]o = const, wi = 0, 



and 



ni (r, (^) = ni[i] (r) + ^ipj (r) J r/3 (r, v?) / [^\JVhM^\) ,^4^0; 

= ^i[i] (^) + ^i[2] {r) j dip ri3 (r, (fi),ril = 0; (2) 
= ni[i] (r) + r2i[2] (r) ^ d^/ (y\r]JjVf)Tj ,^3 = 0; 

where the limit d^\r]4\/dip — >■ is considered for e ^ and the functions 773 (r, ip) and 
'^i[i,2] are stated by some boundary conditions, there is defined a static black hole 
with ellipsoidal horizon (a black ellipsoid). 

Another class of black ellipsoids with anisotropic conformal symmetries, with a 
nontrivial conformal factor Vt which induces a nonzero value for the coefficient wi (see 
details in Ref. [S]), can be defined by a metric (Q) with the data 

VL{r,Lp) = 1 — er'^ sin^ 9ri3{r, 9, ip); 

ri{r,ip) = r'r]^(r,9,(p)sm^ 6 (r"^ ~ 2mr^ /2m; 

ewi{r,(p) = diil{d^l/d(p)~^,dri3/d(p ^ 0, 

= 0,dr];/dip = 0; 

and 

ni(r, (p) = (r) + riip] (r) / dip 773 (r, ip) , (3) 



where, for this type of solutions, the function rj^ (r, 9, ip) is chosen as to have a value 
0(r, ip) = 1]^ (r, 9, ip) sin^ 9 depending only on variables r and ip. 

The maximal analytic extensions of such locally anisotropic black hole metrics in 
the framework of general relativity theory with anholonomic frames were constructed 
in Refs. 12 El- It was shown that the condition of vanishing of the coefficient 1 — 
2m/r + £r){r, ip)/r'^ before 5t defines the equation of a static non-spherical horizon for 
a so called "locally anisotropic" black hole (for a corresponding parametrization of t] we 
may construct a resolution ellipsoid horizon). There is a similarity of the metrics of type 
(QJ with the Reissner-Norstrom solution: the anisotropic metric is obtained as a linear 
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approximation on e from some exact vacuum solutions of the Einstein equations, but 
the Reissner-Norstrom one is a static exact solution of the Einstein-Maxwell equations. 

We can apply the perturbation theory for the metric ((H) (not paying a special 
attention to some particular parametrization of coefficients for one or another class 
of anisotropic black hole solutions) and analyze its stability by using the results of 
Ref. j4j for a fixed anisotropic direction, i. e. by imposing certain anholonomic frame 
constraints for an angle ip = LfQ but considering possible perturbations depending on 
three variables {u^ = = r, = = 9, = t). If we prove that there is a 
stability on perturbations for a value (fo, we can analyze in a similar manner another 
values of (p. A more general perturbative theory with variable anisotropy on coordinate 
if, i. e. with dynamical anholonomic constraints, connects the approach with a two 
dimensional inverse problem which makes the analysis more sophisticate. 

We note that in a study of perturbations of any spherically symmetric system and, 
for instance, of small ellipsoid deformations, without any loss of generality, we can 
restrict our considerations to axisymmetric modes of perturbations. Non-axisymmetric 
modes of perturbations with an e*""^ dependence on the azimutal angle ip {n being an 
integer number) can be deduced from modes of axisymmetric perturbations with n = 
by suitable rotations since there are not preferred axes in a spherically symmetric 
background. The ellipsoid like deformations may be included into the formalism as 
some low frequency and constrained petrurbations. 

For simplicity, in this paper, we restrict our study only to fixed values of the coordi- 
nate (p assuming that anholonomic deformations are proportional to a small parameter 
e; we shall investigate the stability of solutions only by applying the one dimensional 
inverse methods. The metric (pQ) to be investigated is with a deformed horizon and 
transforms into the usual Schwarszschild solution at long radial distances (see Refs. 
[21 El for details on such off-diagonal metrics). 

Let us consider a quadratic metric element 



5u^ 




(4) 



where 




(5) 



qi{x\t) 



with 




(6) 



fi(r, v9o)r■^ 
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and some non-trivial values for qf^ and enj, 



ni = e ni[i](r) +ni[2](r) / r]3{r,ip)dp) . 







We are distinguishing two types of small deformations from the spherical symmetry. 
The first type of deformations, labeled with the index (e) are generated by some e- 
terms which define a fixed ellipsoid like configuration and the second type ones, labeled 
with the index (^), are some small linear fluctuations of the metric coefficients 

The general formluas for the Ricci and Einstein tensors for metric elements of class 
Q with rii = are given in We compute similar values with respect to anholnomic 
frames, when, for a conventional splitting u°' = {x\y"'), the coordinates and y°' 
are treated respectively as holonomic and anholonomic ones. In this case the partial 
derivatives d/dx^ must be changed into certain 'elongated' ones 
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dx^ 


op 





see details in Refs [H 12 El- In the ansatz (jH), the anholonomic contributions of Wi 
are included in the coefficients qi{x'',t), there is only one nonzero value Wi and, in 
consequence, we have to introduce elongations of partial space derivatives only on the 
t-variable. For convenience, in the Appendix we present the necessary formulas for 
Rap (the Ricci tensor) and Ga/s (the Einstein tensor) computed for the ansatz (jlj) with 
three holonomic coordinates (r, 6, p) and on anholonomic coodinate t (in our case, 
being time like), with the partial derivative operators 

Odd d 
di^ Si = — - ni — ,d2 = ^,83 = 

or Ot Od Op 

and for a fixed value p^. 

A general perturbation of an anisotropic black-hole described by a quadratic line 
element (jH) results in some small quantities of the first order uj and g^, inducing 
a dragging of frames and imparting rotations, and in some functions /Xq, with small 
increments 5^ai which do not impart rotations. Some coefficients contained in such 
values are proportional to e, another ones are considered only as small quantities. The 
perturbations of metric are of two generic types: axial and polar one. We investigate 
them separately in the next two Sections. 



3 Axial metric perturbations 



Axial perturbations are characterized by non-vanishing uj and qi which satisfy the 
equations 

R^i = 0, 
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see the explicit formulas for such coefficients of the Ricci tensor in the Appendix. The 
resulting equations governing axial perturbations, SR31 = 0, SR32 = 0, are respectively 



6, (M'^^'-^'-^'Qu 



(7) 



where 

Qij = Siqj - 6jqi, = d^qi - 6iU (8) 

and for /Xj there are considered unperturbed values fifK Introducing the values of 
coefficients (jSj) and (jHl) and assuming that the perturbations have a time dependence 
of type exp(zo"t) for a real constant cr, we rewrite the equations ((H) 



1 + e (A^i + 3r20/2) 



A 



sin^ Or]^^"^ 



5i\Q 



sin^ Orj^^'^ 
T , ^ 1 2 



—ia5rUJ — a'^qi, (9) 
iadeuj + a'^q2 (10) 



for 

Q^") (r, e, (^0, t) = AQ12 sin^ 6 = A sin^ ^(Sagi - 5ig2), A = - 2mr, 

where cj) = for solutions with i7 = 1 and 0(r, yj) = 773 (r, 6', yj) sin^ ^, i. e. r^a (r, 6*, ip) ~ 
sin~^ ^ for solutions with fl = 1 + e.... 

We can exclude the function to and define an equation for Q^"^^ if we take the sum of 
the Q subjected by the action of operator 82 and of the (fTUj) subjected by the action 
of operator 61. Using the relations (jH)), we write 



r% 



sin^ 682 



A 3/2 



61 



+ 



1 + e(A-i + 3r20/2) 

• 3 /) 3/2 
sm t^?73 



+ 



^2^,4 



Q 



0. 



The solution of this equation is searched in the form Q^"^^ = Q + eQ^^^ which results in 



r'^81 



A 



4 3/2 



8,Q]+ sin^ 682 



■ 3 /) 3/2 
,sm urj^ 



d2Q 



Q = eA (r, 9, ^0) 



'111 



where 



- sm-" 682 



l + e(A-i + 3rV/2) 



• 3/33/2 
sm Urj^ 



52Q« - 



with a time dependence like exp[i(Tt] 
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It is possible to construct different classes of solutions of the equation (fTT|) . At the 
first step we find the solution for Q when e = 0. Then, for a known value of Q (r, 9, (po) 
from 

Q('') = g + £Q«, 

we can define Q*-^-* from the equations (jH)) and (fTO|) by considering the values propor- 
tional to e which can be written 



Bi (r, 6*, v3o) , 
B2 (r, e, ¥?o) , 



(12) 



where 



and 



Bi = ni 



dQ 
dt 



Q 



A 



— r 



B2 = -{A-^ + 3r^^/2)d2Q - aV^ sin^ 9 wi. 

The integrability condition of the system ()12|) . 9ii?2 = c^2-Bi imposes a relation between 
the polarization functions r]3,ri,wia.nd rii (for a corresponding class of solutions, see 
formulas (j2I), or (jS))). In order to prove that there are stable anisotropic configurations 
of anisotropic black hole solutions, we may consider a set of polarization functions when 
A (r, 9, (po) = and the solution with Q*-^-* = is admitted. This holds, for example, if 

Aril = "nQr^rjl^"^, Uq = const. 
In this case the axial perturbations are described by the equation 



3/2 4o 
V3 r 



4 3/2 



diQ + sin^ 9d2 



f 1 



Vsin^ 9 



d2Q 



A 



-Q = 



(13) 



which is obtained from (fTT|) for -q^ = rj^ (r, po) , or for 0(r, Pq) = 773 (r, 9, po) sin^ 9. 

In the limit rj^ ^ 1 the solution of equation ()13|) is investigated in details in Ref . jl] . 
Here, we prove that in a similar manner we can define exact solutions for non-trivial 
values of rj^. The variables r and 9 can be separated if we substitute 

Qir,9,po) = Qo{r,Po)C[^TiO), 

where are the Gegenbauer functions generated by the equation 



^ sin^^ 9-^ + n{n + 2u) sm^" 9 
d9 d9 



C:{9) = 0. 



The function Ci_^2'^{9) is related to the second derivative of the Legendre function Pi{9) 
by formulas 



• 3 /I 



1 dPi{9) 
sin 9 d9 



The separated part of ()13|) depending on radial variable with a fixed value ip^ transforms 
into the equation 

li'^A' 



A ^ 



A dQo 



Qo = o, 



(14) 



where /i^ = (/ — + 2) for / = 2,3, ... A further simphfication is possible for 773 
773 (r, lpq) if we introduce in the equation (jl4j) a new radial coordinate 



and a new unknown function Z^^^ = r ^Qo('")- The equation for Z'^'^'^ is an Schrodinger 
like one-dimensional wave equation 



+ ) Z^'') = (15) 



rlr"^ 3/2 



with the potential 



A 



5 3/2 



2 4 / A 
/X — r — 



(16) 



and polarized parameter 



5' = aVr]: 
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This equation transforms into the so-called Regge-Wheeler equation if 773 = 1. For 
instance, for the Schwarzschild black hole such solutions are investigated and tabulated 
for different values of Z = 2, 3 and 4 in Ref. 

We note that for static anisotropic black holes with nontrivial anisotropic conformal 
factor, Q = 1 + e..., even 773 may depend on angular variable 6 because of condition that 
0(r, ipo) = rjs (r, 9, ipo) sin^ 9 the equation (fT^ transforms directly in (fT3j) with /i = 
without any separation of variables r and 9. It is not necessary in this case to consider 
the Gegenbauer functions because Qo does not depend on 9 which corresponds to a 
solution with / = 1. 

We may transform ()15|) into the usual form, 

^ + aA Z^") = V^^)Z'^'^^ 



dr. 
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if we introduce the variable 



J c/r#r/3^/^(r#,(^o) 



for V^"^^ = Tj^^'^V'^'^^ So, the polarization function 773, describing static anholonomic 
deformations of the Scharzshild black hole, " renormalizes" the potential in the one- 
dimensional Schrodinger wave-equation governing axial perturbations of such objects. 

We conclude that small static "ellipsoid" like deformations and polarizations of 
constants of spherical black holes (the anisotropic configurations being described by 
generic vacuum off-diagonal metric ansatz) do not change the type of equations for 
axial petrubations: one modifies the potential barrier, 

\/(-) = ^ [(^/i2 ^ 2) r - 6m] — > V^"^^ 

and re-defines the radial variables 

r^, = r + 2m In (r / 277i — 1) — > ((/Pq) 

with a parametric dependence on anisotropic angular coordinate which is caused by 
the existence of a deformed static horizon. 
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4 Polar metric perturbations 

The polar perturbations are described by non-trivial increments of the diagonal metric 
coefficients, = 5fj,^^^ + Sfj,^^^ , for 

where 5jj!£>{x'^ ^t) parametrize time depending fluctuations which are stated to be the 
same both for spherical and/or spheroid conflgurations and is a static deformation 
from the spherical symmetry. Following notations (0) and (jH)) we write 



and 



= r/sJ\Kle''^ = r^J\q{r)\,e''' = r/ig sin^, e^* = A/r^ 



where = for the solutions with Q = 1. 

Examining the expressions for i?4j, R12R33 and Gu (see the Appendix) we conclude 
that the values Qtj appear quadratically which can be ignored in a linear perturbation 
theory. Thus the equations for the axial and the polar perturbations decouple. Con- 
sidering only linearized expressions, both for static e-terms and fluctuations depending 
on time about the Schwarzschild values we obtain the equations 



61 (5/i2 + 5/i3) + (r-i - 5i/i4) (5/i2 + 5/i3) - 2r~^6i2i = (5i?4i = 0) , 
82 {6fii + 5/i3) + ((5/i2 - ^/is) cot 6' = {6R42 = 0) , 

£^2^1 (5/^3 + Sfli) - 5i (5/i2 - 5/i3) cot 6- 

(r^^ - (5i/i4) d2{6pii) - [r^^ + 61^^ d2{6^i) = {6R42 = 0) , 
e2/^4|2 (^-1 + s,fi,) 6,{6fi3) + r~'6, (5/^3 + 6fi, - 6fi, + 6(12) + (17) 
Si [SiiSfis)] - 2r-'6fi, (r-i + 251^4)} - 2e-^^'d,[d4{S^i3)] + 

r-^{d2[d2{6fi3)] + d2{26fi3 + SiJ4 + Si2i~6fi2)cote + 26iJ2} = (5i?33 = 0) , 

2r-'6fii (r-i + 26ifii)] - e-^^^d^[d^iSfi3 + 5/^2)] 
+r-^{d2[d2{Sfi3)]+d2{26fi3 + Sfi4-Sfi2)cote + 26fi2} = {6Gu = 0) . 

The values of type SjJa = Sfx^^^ + 5jj!£> from (fT7|) contain two components: the first 
ones are static, proportional to £, and the second ones may depend on time coordinate 
t. We shall assume that the perturbations 811^^^ have a time-dependence exp[(T)f:] so 
that the partial time derivative "94" is replaced by the factor ia. In order to treat both 
type of increments in a similar fashion we may consider that the values labeled with 
{e) also oscillate in time like exp[cr*-'^''t] but with a very small (almost zero) frequency 
cr(^) 0. There are also actions of "elongated" partial derivative operators like 
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To avoid a calculus with complex values we associate the terms proportional eriidi to 
amplitudes of type eiriidi and write this operator as 



For the "non-perturbed" Schwarzschild values, which are static, the operator 6i reduces 
to di, i.e. 6iVa = diVa- Hereafter we shall consider that the solution of the system (fTTjl 
consists from a superposition of two linear solutions, 6^a = ^/^a^ + ^IJ-a^', the first 
class of solutions for increments will be provided with index (e), corresponding to 
the frequence a^^^ and the second class will be for the increments with index (<^) and 
correspond to the frequence a^''\ We shall write this as S/j^^ and a^A) for the labels 
A = e or q and suppress the factors exp[(T*^^^t] in our subsequent considerations. The 
system of equations (jl7|) will be considered for both type of increments. 

We can separate the variables by substitutions (see the method in Refs. [311]) 



T^^) (r)Pi (cos + V^^\r) cot 9dPi/d9 
and reduce the system of equations (fT7|) to 



(19) 



and 



2r'Hi (iV(^)) - /(/ + l)r-2e-2^*iV(^) - 2r-\r-^ + 2diUi)L^^^ - 2{r-'+ (20) 



diu^)6i [n^^^ + (/ - 1)(/ + 2)V(^V2] - (/ - 1)(/ + 2)r-^e^^'"' (v^^^ - L^^^) - 

L(^) + (/- l)(/ + 2)\/(^V2 



where we have introduced new functions 

X(^) = ^(/-l)(/ + 2)\/(^) 

and considered the relation 

T(^) - r(^) + l(^) = {6R,2 = 
We can introduce the functions 



0, 



0). 



L^^^+ea^A) 



niL^^Ur, Ar(^) = A^^^) + ea^A) [ riiN^^Ur, (21) 



NiA)^ y{A) ^ y(A) ^ 



n 



LiA) 

f{A) 

for which 
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and, this way it is possible to substitute in (fT^ and (PU]) the elongated partial derivative 
6i by the usual one acting on "tilded" radial increments. 

By straightforward calculations (see details in Ref. 0]) one can check that the 
functions 



,6mX(^)/r(/ -!)(/ + 2) -L(^) 



' r(/- l)(/ + 2)/2 + 3m 
satisfy one-dimensional wave equations similar to (fT3j) for Z^^^ with 773 



1, when 



+ o-f,, 1 



(+) 

(A) j ^{A) 
y{+) 



T/(+)v(+) 



(22) 



where 



2A 



r5[r(/- l)(/ + 2)/2 + 3m]^ 



X {9m^ 



2^ 



+ l(/-lf(/ + 2)V 



l + -(/-l)(/ + 2) + 



l)(/ + 2) + m 
3m1 



(23) 



For e — > 0, the equation transforms in the usual Zerilli equation [TT| 1^. 

To complete the solution we give the formulas for the "tilded" X- and A^- 
f actors, 



1(A) 



(/-!)(/ + 2) 



($ 



i)([+2)~(+) 
2r 

(A) 



(24) 



2r 



m 



(A) 



(/-l)(/ + 2)r 



r — 2m y 
(/-!)(/ + 2) 



2(/ -!)(/ + 2) + 12m dr# 



[r{l- 
12m2 



l)(/ + 2) + 6m]^ 



X 



3m(/ -!)(/ + 2) + -(/ 



l)(/ + 2) [(/-!)(/ + 2) + 2] 



where 



$(^) = (/_i)(/ + 2)e"^ y 



?(+) 



-dr. 



(/-l)(/ + 2)r + 6m 

Following the relations ()21|) we can compute the corresponding "untileded" values an 
put them in (|T8|) in oder to find the increments of fluctuations driven by the system of 
equations ()17j) . For simplicity, we omit the rather compersome final expressions. 

The formulas ()24j) together with a solution of the wave equation (j^^ complete the 
procedure of definition of formal solutions for polar perturbations. In Ref. |4j there are 
tabulated the data for the potential (j^Hj) for different values of / and (/ — 1) (/ + 2)/2. In 
the anisotropic case the explicit form of solutions is deformed by terms proportional to 
eriia. The static ellipsoidal like deformations can be modeled by the formulas obtained 
in the limit (T(e) — > 0. 
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5 The Stability of Black Ellipsoids 



The problem of stability of anholonomically deformed Schwarzschild metrics to external 
perturbation is very important to be solved in order to understand if such static black 
ellipsoid like objects may exist in general relativity. In this context we address the 
question: Let be given any initial values for a static locally anisotropic configuration 
confined to a finite interval of r^, for axial perturbations, and r^, for polar peturbations, 
will one remain bounded such peturbations at all times of evolution? 

We have proved that even for anisotropic configurations every type of perturbations 
are governed by one dimensional wave equations of the form 

3^ Z o 

—— + a^Z = VZ 25 
dp 

where p is a radial type coordinate, Z is a corresponding Z*^^) or zjj^j with respective 

smooth real, independent of cr > potentials V^^^ or V^~^ with bounded integrals. For 
such equations a solution Z{p, cr, ipo) satisfying the boundary conditions 

T{a)e'''P (p ^ -oo) 

(the first expression corresponds to an incident wave of unit amplitude from +oo giving 
rise to a reflected wave of amplitude -R(cr) at +oo and the second expression is for a 
transmitted wave of amplitude T(cr) at — oo), provides a basic complete set of wave 
functions which allows to obtain a stable evolution. For any initial perturbation that 
is smooth and confined to finite interval of p, we can introduce the integral 

^(p,0) = (27r)^^/' r^^{a,0)Z{p,a)da 

J — oo 

and define, at later times, the evoluiton of perturbations, 

ij{p, t) = (27r)"^/^ / ij{(T, 0)e*"*Z(p, (T)d(T. 



The Schrodinger theory garantees the conditions 

mp,0)\'dp= / \ij{a,0)\'da= / mp,0)\'dp, 

-oo J —oo J —oo 

from which the boundedness of ip{p,t) follows for all t > 0. 

In our consideration we have replaced the time partial derivative d/dt by io", which 
was represented by the approximation of perturbations to be periodic like e^"^. This is 
connected with a time-depending variant of (j^ . like 

df^ dp^ 
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Multiplying this equation on dZjdt, where Z denotes the complex conjugation, and 
integrating on parts, we obtain 

+^(dZd'^Z dZ d'^Z ^^BZ\^ 

providing the conditions of convergence of necessary integrals. This equation added to 
its complex conjugate results in a constant energy integral. 



dZ 



dt 



+ 



dZ 



dp 



+ V\Z\ \ dp = const, 



which bounds the expression \dZ/dt\'^ and excludes an exponential growth of any 
bounded solution of the equation (j2Sl)- We note that this property holds for every 
type of "ellipsoidal" like deformation of the potential, V ^ V + eV^^\ with possible 
dependences on polarization functions as we considered in (fT^ and/or ((221) • 

The general properties of the one-dimensional Schrodinger equations related to per- 
turbations of holonomic and anholonomic solutions of the vacuum Einstein equations 
allow us to conclude that there are locally anisotropic static configuratios which are 
stable under linear deformations. 

In a similar manner we may analyze perturbations (axial or polar) governed by a 
two-dimensional Schrodinger wayve equation like 

d^Z d^Z ^, ^d^Z 

+ A{p,^,t)—r-V{p,^,t)Z 



dt"^ dp^ 



for some functions of necessary smooth class. The stability in this case is proven if 
exists an (energy) integral 



TT r + oo 



dZ 



dt 



dZ 



dp 



+ 



^ dp 



+ V\Z\ dpd(p = const 



which bounds {dZ/dtl"^ for two-dimensional perturbations. For simplicity, we omitted 
such calculus in this work. 

Finally, we note that this way we can also prove the stability of perturbations along 
"anisotropic" directions of arbitrary anholonomic deformations of the Schwarzschild so- 
lution which have non-spherical horizons and can be covered by a set of finite regions 
approximated as small, ellipsoid like, deformations of some spherical hypersurfaces. 
We may analyze the geodesic congruence on every deformed sub-region of necessary 
smoothly class and proof the stability as we have done for the resolution ellipsoid hori- 
zons. In general, we may consider horizons of with non-trivial topology, like vacuum 
black tori, or higher genus anisotropic configurations. This is not prohibited by the 
principles of topological censorship if we are dealing with off-diagonal metrics and 
associated anholonomic frames The vacuum anholonomy in such cases may be 
treated as an effective matter which change the conditions of topological theorems. 
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6 Outlook and Conclusions 



It is a remarkable fact that, in spite of appearance complexity, the perturbations of 
static off-diagonal vacuum gravitational configurations are governed by symilar types 
of equations as for diagonal holonomic solutions. The origin of this mystery is located 
in the fact that by anholnomic transforms we effectively diagonalized the off-diagonal 
metrics by "elongating" some partial derivatives. This way the type of equations gov- 
erning the perturbations is preserved but, for small deformations, the systems of linear 
equations for fluctuations became "slightly" nondiagonal and with certain tetradic 
modifications of partial derivatives and differentials. In details, the question of relating 
of particular integrals of such systems associated with systems of linear differential 
equations is investigated in Ref. In our case one holds the same relations between 
the potentials K*-^^ and V^*^~) and wave functions Z^^^ and Z^^-j* with that difference 
that the physical values and formulas where polarized by some anisotropy functions 
ri3{r,9,{p),Q{r,{p),q{r),ri{r,ip), Wi{r,{p) and ni(r, and deformed on a small param- 
eter e. 

We also observe that the "anisotropic" potentials V^^^ and V^~^ may be defined for 
such polarizations ri'i{r,9,ip) as they would be smooth functions, integrable over the 
range of r^,, (—00, +00) and positive everywhere. For real cr, the anisotropic solutions, 
represent ingoing and outgoing waves of type e^*°"''*(r^, — > ±00). We conclude that 
the underlying physical problem is one of reflexion and transmission of incident waves 
(from -|- or — 00) by some anisotropically deformed one-dimensional potential barriers 
V^'') and/or V^~\ for every fixed anisotropic angle <^o- Because at distancies fare away 
from deformed horizons, the metrics of anisotropic black holes transforms into the usual 
Schwarzschild solution we can prove the equality of the reflexion and the transmission 
coefficients for the axial and the polar perturgations, even the wave equations are an- 
holonomically deformed. Such formulas are tabulated in Ref. ^ and may be used 
for anisotropic solutions but with some redefined coefficients. The general properties 
of the one-dimensional potential-scattering are preserved for small anholonomic de- 
formations; they are concerned with some solutions of one dimensional Schrodinger's 
wave equations. 

We emphasize that we proved the stability of black ellipsoids by fixing any anisotro- 
pic directions in the off-diagonal metrics, if = ip^ i. e. under certain anholonomic 
constraints imposed on vacuum gravitational configurations. A more general con- 
sideration with variable (p relates the problem to the two dimensional Schrodinger's 
potential-scattering problem, as well to an anholonomic Newman-Penrose formalism 
which makes the solution of the problem of stability to be more sophisticate. Neverthe- 
less, we are having strong arguments that the stability of anisotropic static solutions 
proved in the parametric approximation ip = ipQ will hold at least for such polarizations 
which model two dimensional scattering effects containing some particular one dimen- 
sional stable Schrodinger's potentials. The one-dimensional perturbation analysis is 
the first, very important, step in investigating the stability of any physical system; it 
should be included into the higher dimensional and/or less constrained approaches. 

We found that the origin of the results that the functions Z^~^ and Z^~^^ and their 
anisotropic extensions Z^'^'^ and zj^} do in fact satisfy very similar one-dimensional 
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wave equations steems from the fact that some deep properties of static solutions of 
the Einstein equations which hold both for diagonal and off-diagonal static-stationary 
vacuum and/or electrovacuum metrics. 

Finally, we conclude that there are static black ellipsoid vacuum configurations 
which are stable with respect to one dimensional perturbations, axial and/or polar ones, 
governed by solutions of the corresponding one-dimensional Schrodinger equations. 
The problem of stability of such objects with respect to two, or three, dimensional 
perturbations, and the possibility of modeling such perturbations in the framework 
of a two-, or three-, dimensional inverse scattering problem is a topic of our further 
investigations. 
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Appendix 

We compute the coefficients of the curvature tensor as 

T)Ct e -pa e -pa . -pa per pa per 

of the Ricci tensor as 

and of the Einstein tensor as 

Gj3-f = -R/37 — -gp'^R 

for R = gf^'^Rp^ and 5^ = d-y for 7 = 2,3,4 and 61 = di — riid/dt. Straightforward 
computations for the quadratic line element (jH) give 

^11 = -e"^'''[5ii(/i3 + Ai4 + Ai2) + ^lyUs^il/^s - Ail) + 5iAi25i(Ai2 - Ail) + 
5i/i4(5i(Ai4 - A^i)] - e~^^^[d22fii + d2fiid2{fi3 + At4 + A^i - At2)] + 

e-^'^'ldlfii + d,fiid^{fi3 - /i4 + A^i + Ai2)] - ^e^^'^'-^'^e-^'^'Ql^ + e 

R12 = -e~^^'^''[d2Si{fi3 + fi2) - d2filSi{fl3 + fii) - 6ifi2d4{fi3 + ^J'l) 
+Slfl3d2fi3 + 5l/i452Ai4] + ^e^^'-^'''-'''-'''Ql4Q24, 

^33 = -e"^''' [SufJ'S + SiHsSi (fis + At4 + At2 - A^i)] - 

e~'^'[di2f^3 + d2l^3d2{fi3 + AX4 - AX2 + Ail)] + ie2(^«-^i-^2)Q2^ + 
e-'^'[dlfi3 + d^fisd^ifis -fi4 + fi2 + A^i)] - ^e^^^'->^'\e-'^'Ql, + e 
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Ril = -e '"'[di5i{fX3 + 1x2) + SifXsdiiiJ^s - fJ'i) + SiiJ.2di{iJ.2 - fJ'i) 
-R44 = -e"^^^ [d^ifi'i + /i2 + yWs) + d^iJ^sd^if^s - ^14) + difiid^ifii - fi^) + 

e-^^'[dl2li, + a2/i492(/i3 + /i4 - /il + /i2)] - ^e2('^3-/.4) [^-2/.! g2^ ^ C-^^^QI,], 

where the rest of coefficients are defined by similar formulas with a corresponding 
changings of indices and partial derivative operators, i?22, -R42 and R32 is like R^i 
and i?3i with with changing the index 1 — > 2. The values Qij and (5j4 are defined 
respectively 

Qij = Sjqi - 6iqj and Qi4 = d^qt - 6iUJ. 
The nontrivial coefficients of the Einstein tensor are 

Gu = e"^'"" [<922(At3 + /i4) + C^2(/i3 + /i4)c^2(/i4 " /i2) + ^^2/^3^2/^3] " 
e~^''*[<944(/i3 + H2) + d^ifXs + /i2)54(/i2 - /X4) + 54/^3<94Ai3] + 
e"^''^ [5l/i4 + 5l (/U3 + yU2) + 5l/i3^lAt2] - 

lg2AJ3 |-g-2(/ii+fi2)Q2^ — e~^*-^^^''''^(5i4 + e~^'-'^^''''^^^Q24] 

^33 = e~^''45ii(/i4 + /i2)+5l/i45l(/i4-yUl+/i2)+5l/i25l(/i2-/Ul)] + 
e"^'''[522(/U4 + /il) + 92(1X4 - /i2 + /il) + 52/il52(/Ul " /U2)] " 
e~^''*[(944(/ii + /i2) + dijJidiifii - /i4) + (94/^2C^4(yU2 - /i4) + di^lidi^l2] + 

-p2M3rp-2(m+M2)r)2 _ p-2{m+M4)r)2 _ p-2(M2+A«3)n2 i 
[e i<^/]^2 Vi4 V24J5 

G44 = e"2^i[5^^(/i3 + ;U2) + (5i/i35i(/i3 - ;Ul + /i2) + (5i/i25l(Ai2 - /il)] - 

e~"^n522(/"3 + /^i) + ^2{^i3 - /X2 + yui) + 52/iia2(/ii - ^^2)] - ^e^^'^^-'^^-^^^Q^ 
+e-2'^^[a4/.3a4(^i + /i2) + 54^i94/i2] - -^6^^^^^- ^^\e-^^^^ Q\, - e~^^^Ql^. 
The component G22 is to be found from Gu by changing the index 1^2. 
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